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Abstract 

A dynamical model is proposed to describe tlie coupled decomposition and profile evolution of a free 
surface film of a binary mixture. An example is a thin film of a polymer blend on a solid substrate under- 
going simultaneous phase separation and dewetting. The model is based on model-H describing the cou- 
pled transport of the mass of one component (convective Cahn-Hilliard equation) and momentum (Navier- 
Stokes-Korteweg equations) supplemented by appropriate boundary conditions at the solid substrate and the 
free surface. 

General transport equations are derived using phenomenological non-equilibrium thermodynamics for 
a general non-isothermal setting taking into account Soret and Dufour effects and interfacial viscosity for 
the internal diffuse interface between the two components. Focusing on an isothermal setting the resulting 
model is compared to literature results and its base states corresponding to homogeneous or vertically 
stratified fiat layers are analysed. 
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I. INTRODUCTION 



Driven by applications in coating technology, micro and nano structuring of soft matter layers, 
and, in general, the development of (multi-)functional surfaces the understanding of thin films of 
simple and complex fluids is of growing importance. Recent years have seen on the one hand 
major advances in experimental techniques of preparation and analysis and on the other hand 
intense developments of the theoretical description of the statics and dynamics of homogeneous 
and structured fiinis.-Ja2i2ii^i^il2i2 The dynamics of the structuring often represents examples for 
micro- and nano-fluidic flows, a present focus of interest in its own right. 

For thin one-layer free surface films of one-component simple or polymeric liquids experi- 
mental results2ii^»^^ and theoretical understandin g^'*' ^^I'^i ' 8, 19,20 ^^^l developed. However, 
experiments increasingly focus on complex situations like the evolution of multilayer films of 
partially miscible^' or immiscible"^ liquids, complex fluids like polymer blends that might undergo 
dewetting or/anJ decomposition^^-^^i^l or solutions of polymers, nanoparticles, colloids or polymer 
blends with interacting convective motion, phase separation, evaporation/condensation and evolv- 
ing rheology .^'*i^^i^^'^"^i^^i^^i^°'^'i^^i^^i^"^'^^i^^i-^"^i^^'^^ Theoretical descriptions exist, however, only for a 
small part of the experimentally known complex scenario and phenomena involving free surface 
thin films. Recent advances include a fully nonlinear thin film description in long wave approxima- 
tion for two layers of immiscible liquids under air^^iiLl^idSjii and between two plates,— the analysis 
of the dewetting behaviour on chemically or topographically heterogeneous s ub strates ^46,47.48,49 j-j^g 
study of the dynamics of depinning of a driven drop on a heterogeneous substrate,— the description 
of films with surface active nanoparticles,— the inclusion of evaporation/condensation in the thin 
film description .'^i^^i^^i^^'^'^i^^i^^ 

Thin films of polymer blends are one of the 'simplest' complex systems listed above and 
extensive experimental results can be found in the literature .'^i^^i^^i^"^i^^i^^i^°i^^i^^i^^i^'^-^^'^^i^"^ De- 
scribed effects include the dependency of the evolving structure (stratified bi- or tri-layer struc- 
ture that is laterally homogeneous; purely lateral phase separation, checkerboard structure) 
on substrate properties, -^^^^^Mi^ surface roughening or film morphology changes during phase 
separation ,^^i^^i^^i^° surface directed spinodal decomposition,— subsequent vertical phase separa- 
tion and dewetting,-^i^ and surface phase inversion.— The influence of heterogeneous substrates 
was also studied.'^^ However, as detailed below to our knowledge there exists no theoretical de- 
scription of the involved processes that takes into account the evolving composition of the mixture 



2 



and the evolving surface profile of the film. 

The aim of the present paper is to present such a description based on the underlying transport 
equations. To cover the coupled time evolution of the film thickness and concentration profiles one 
has to supplement the coupled transport equations for momentum and concentration by appropriate 
boundary conditions at the free surface and at the solid substrate. 

Two groups of studies are present in the literature that address part of the involved questions. 
On the one hand, the classical Cahn-Hilliard model^^ describing purely diffusive decomposition 
of a binary mixture was studied for films in a gap between two solid plates .-^^i^^ At the plates 
boundary conditions prescribe zero diffusive flow through the plates, energetic preference of the 
plates for one component and an enforced (or reduced) demixing at the plates. However, such 
a model can in principle not account for an evolving surface deflection in decomposing films as 
observed in phase-separating polymer blends on homogeneous^^*^^*^^ or patterned^^ substrates. 
An evolving film profile is by definition related to a convective flow of the mixture. 

On the other hand, the coupling of momentum and concentration transport for weakly miscible 
fluids (i.e. decomposing mixtures) is well studied for bulk systems using the so-called model- 
76,77,78,79,80 jj- couplcs the convcctivc Cahn-Hilliard equation and the Navier-Stokes equations 
ammended by a concentration dependent stress tensor, the so called Korteweg stresses.*^' For a 
survey of the history see Ref.S. A variant of model-H is also used for the dynamics of momentum 
and density of a single component near a liquid-gas phase transition in isothermal^^'^^iMiS^ or non- 
isothermal conditions .■^^'^^ Two-phase liquids and binary mixtures between two solid plates are 



also investigated . ^ ^ A similar model for miscible liquids was studied, for instance, in Refs. 
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and|90l There the convective Cahn-Hilliard equation is replaced by a 'normal' convective diffusion 
equation. Nowadays, model-H is also applied to multiphase flows in closed micro channels of 
different geometries (straight quadratic channel,'^! T-junction.—) 

To describe a film of a mixture under air, model-H is used to describe the dynamics within 
the film, i.e. modelling the creation and evolution of 'internal' diffuse interfaces within the film. 
The bulk model has to be completed by boundary conditions at the solid substrate (discussed in 



Refs. 



91 



94 



1951 ) and by boundary conditions at the free surface. The latter represents an 'external' 
sharp interface. The interaction of the internal diffuse interfaces with the external sharp one via a 
solutal Marangoni effect results in an additional driving force of the evolution. 

Other related work involves an ad-hoc lubrication approximation model coupling evolution 
equations for the film thickness and the mean concentration in the film.'^^ We believe that such 
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a model might be correct in the limit of weak vertical variation of concentration but is not able to 
describe vertically stratified films and their evolution. An alternative approach uses microscopic 



discrete models like lattice gas (for results on a binary alloy see Ref. |98|) or Molecular Dynamics.— 
This leads, however, often to a strong restriction in length- and time-scales that can be studied. 

We present our work as a sequence of papers. The present first part derives and discusses 
the basic transport equations and analyses steady base states. The accompanying second part 
performs a detailed stability analysis with respect to transversal instability modes for the various 
qualitatively different base states. Thereby the consequences of convective transport are studied in 
detail and the sequence of patterning processes is predicted. A planned further sequel will focus 
on the nonlinear evolution. 

The present paper is structured as follows. In Section |ll] we derive the coupled transport 
equations for momentum, density and temperature in the framework of phenomenological non- 
equilibrium thermodynamics. After discussing the physical interpretation of the individual contri- 
butions representing, for instance, an internal Soret effect (i.e. Soret effect for the internal diffuse 
interface) and internal interface viscosity, the model is simplified assuming an isothermal setting, 
vanishing interface viscosity and internal energies resulting from a setting close to the critical point 
of demixing. The resulting model-H is compared to versions found in the literature discussing the 
issue of defining pressure and chemical potential. Section|III]introduces boundary conditions at the 
rigid solid substrate and the free interface. It is explained in detail why the incorporation of con- 
vective flow is a necessary precondition to describe evolving surface deflections. Next, Sections ITVl 
and |V] introduce the non-dimensionalization and local energies, respectively. Homogeneous and 
vertically stratified, transversally homogeneous steady state solutions are analysed in Section IVTl 
The final Section IVIIl summarizes, compares to the literature and gives an outlook on the sequel. 
The Appendix lAl uses variational calculus to independently derive the boundary conditions in the 
static limiting case. 

n. DERIVATION OF EXTENDED MODEL-H 

First we present a derivation of an extended model-H that accounts for all cross couplings 
of the transport equations for momentum, concentration and temperature. This includes Soret- 
and Dufour-effects with nonlinear coefficients and interface viscosity for th e diffuse interface. It 



follows in spirit the derivations of the Navier-Stokes equations given in Refs. 



10' 



I 



101 



The online 
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version shows new terms related to the concentration field in red. 



A. General transport equations 

1. Conserved quantities 

Starting point are the transport, conservation and balance laws for the relevant phenomenolog- 
ical thermodynamic entities. In general, we have for a conserved scalar or vector field a(x, t) the 
transport equation 

|a + V-j'„ = (1) 

where is a general flux density that is a vector or second order tensor. Note that a dotless product 
corresponds always to a tensor (or outer) product, whereas a dot '■' product is an inner product 
(resulting in a tensor of the order n — 2, where n is the order of the respective tensor product). The 
contribution by convective transport with the velocity v is expressed explicitly by j'„ = ja + av 
where is the diffusive flux caused by (several) microscopic mechanisms. 

By definition the total mass density p(x, t) is transported by convection only, i.e. the mass 
density flux is 

jp = pv = g (2) 

corresponding to the momentum density g. The transport equation for the density (continuity 
equation) is 

+ V ■ g = (3) 

The density of the momentum g as well as the density of the total energy e are transported by 
convective and diffusive fluxes, i.e. 

|g + v.j; = o (4) 

|e+V.jl = (5) 

where j'^ is the tensor of the momentum flux density and j'^ is the energy flux density. Note that all 
densities are per volume. 

By explicitly denoting the transport by convection as before, = ev + and = vg + a_, 
where q_ is the usual symbol for the diffusive momentum flux (sometimes also called pressure 
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tensor corresponding to the negative of the stress tensor). We will use underlined symbols to 
denote tensors of 2nd or higher order. Eqs. dS]) and (U) result in 

^e + V-(ev) + V-j, = (6) 

d 

and i^g + V ■ (vg) + \/-a = (7) 
at 

respectively. 

Introducing the material time derivative D/Dt = d/dt + v-V one obtains for the velocity 
field [© and ©] 

P^ + V-o: = (8) 

For a binary mixture of fluids a transport equation for the mass density of one of the components 
has to be added beside the one for the total density p. Choosing pi we have 

d 



at 

i.e. 



pi + v-j;^ = o, (9) 



^Pi + V-(piv) + V-jp, =0. (10) 

The density of the other component is p2 = P — pi- These are all the conserved quantities. The 
conservation of angular momentum is guaranteed by the symmetry properties of the stress tensor 
(see below).— Note that all densities used throughout the paper are volume densities. Next we 
discuss the transport equations for non-conserved quantities. 



2. Non-conserved quantities 

The transport equations for non-conserved quantities contain additional source terms, i.e. for a 
general field a one writes 

^ + v-ji = g„ (11) 

where Qa is a possibly space- and time-dependent source density. Relevant non-conserved quan- 
tities are the densities of the internal energy u and of the entropy s. 

For systems with small gradients of concentration and/or temperature the energies do only 
depend on the local fields. For strong gradients, however, this statement does not hold any more 
and the energy of a system will depend also on field gradients. In the present case we consider 
strong density gradients related to diffuse interfaces between different phases. The underlying 
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assumption is that for a demixing system gradients in pi might be much larger than all other 
gradients. The latter enter the theory as parametric dependencies on space only. 

We define the internal energy u{p, pi, s, t) as the thermodynamic equilibrium value for a local 
fluid element, i.e. it shall not depend on gradients. The conserved total energy e, however, shall 
include gradient terms in pi. The relation between the two is 

e = u + ^v' + ^-{Vprf. (12) 

Note that the unit of ^ is [^] = m^/(kgs^). The energy densities have units [e] = [u] = 
Nm/m^ = kg/ (m s^). The transport equation for the internal energy is 

^ + v-j; = g. (13) 

with 

j; = «v + j,. (14) 

For the irreversible processes in question entropy is not conserved. The transport equation for its 
volume density is 

|^ + V-j', = Q. = | (15) 

with j'^ = sv + the total entropy flux density. We write the source density Qg in the usual form 
R/T where R is the so-called dissipation function and T the temperature.— 



B. Determination of thermodynamic forces 

The flux densities a, j^, ju, j^^, js and source terms R/T, remain to be determined. The 
specific transport equation for the internal energy is obtained from the transport equation of the 
total energy Eq. © using Eqs. Q, ©, and ([HI). It reads 

-^u + uV ■ V + V ■ j„ = {-g: + e(Vpi)(Vpi) - i [|(Vpi)' + piApi] 1} : Vv 

-e(Api)V-j,, (16) 

with 

= je - V ■ - ^(Vpi) [pi(V • V) + V ■ jpj . (17) 

The symbol ':' stands for a double inner product, i.e. a : b = . aijhji. 
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The time evolution of the entropy is deduced using a local form of Gibbs relation for each fluid 
element, i.e. from the assumption that small fluid elements are in thermodynamic equilibrium. 
Gibbs relation for a local fluid element of volume V writes 



dU = TdS -pdV + fii dNi + fi2 dN2 (18) 

where U, S, p, (xi, fi2, Ni, N2 stand for internal energy, entropy, pressure, chemical potentials of 
component 1 and 2, and particle numbers of component 1 and 2. The chemical potentials (with hat) 
are related to particle numbers. The relation (fTSi) is transformed expressing extensive variables by 
the corresponding densities using U = uV , Ni = riiV = pi Na/Mi V , N = Ni + N2, P2 = p — Pi, 
and S = sV, where N is the total particle number, p the density of a mixture. Mi the molar mass 
of component i, and Na is the Avogadro number. One obtains 

du = T ds + P2 dp + Pd dpi + {-u + P2P + /^dPi +Ts - p)dV/V, (19) 

The chemical potential p2 = ^2^0/^2 of component 2 and the difference of the chemical poten- 
tials of components 1 and 2 pd = PiNa/Mi — p2Na/M2 are related to volume densities and have 
units [P2] = [pd] = m^/s^. 

Relation (fT9l) is valid for arbitrary local volume V, i.e. one obtains the local Gibbs relation 

du = T ds + p2 dp + Pd dpi (20) 

and the local Gibbs-Duhem relation 

p = -u + Ts + P2P + PdPi- (21) 

Here, we observe that, within the framework of volume density quantities, p behaves as a thermo- 
dynamic potential that is related with u by the Legendre transform (|2T1) . Furthermore, 

Eq. (I20I) is divided by a small time span dt that is, however, large as compared to typical micro- 
scopic time scales yielding 

du ds dp dpi 

-dI = ^Tt+^'Tt+^'^- ^^^^ 

This relation is valid in all local volume elements that might be convected by the flow, i.e. the 
derivatives d/dt correspond to Lagrangian or material time derivatives denoted above D/Dt. Us- 
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ing equations ©, (flOl) . and (fT6l) one transforms (1231) into the wanted form of Eq. (fTSl) 



5s „ 



^ -^ + e(Vpi)(Vpi)+U-ePlApi - ^(Vpi)^) l[>:(Vv) 



^ 1 +Jpi ■ V ( — ^ 



+ j«-V(;^)+j,,-V(^^^^V^) (24) 



with the pressure p and the flux of internal energy j„ given by (|2T1) and (flTI) . respectively. Here 
it is already possible to see the structure of the dissipative contribution to the pressure or stress 
tensor (the part in {} on the r.h.s.). The reversible entropy transport (cp. Eqs. (flSl) and (|24l) ). i.e. 
the entropy flux 

j; = ^v + i + J^(eApi-pd) (25) 

contains the convective transport, the transport via the heat flux and the transport via the diffusive 
flux of species one. Comparing Eqs. (fT5l) and (|24|) allows to identify the source term for the 
entropy. It is related to irreversible processes and written in terms of the dissipation function 

R = |-^ + e(Vpi)(Vpi)+ ("p-ePiApi - \{Vpif \ l| : (Vv) 



+Tj„ • V ( - ) + T j,, ■ V ( ) . (26) 



We can directly deduce the reversible part of the pressure tensor gi = + g^, because only 
the dissipative part of contributes to the dissipation function, i.e. 

= e(Vpi)(Vpi) +(^p- ^pi^pi - f (Vpi)^) I. (27) 

Note that negative of — pi is known as the capillary or Korteweg stress tensor in the 
literature .^^1^^1^'^^^ The dissipative part is also called viscose pressure tensor or friction ten- 
sor. The dissipation function has the structure R = J2a ' ^a, where the ja and Fq, are general 
thermodynamic fluxes and forces, respectively, that might be tensors. Correspondingly the '■' 
stands here for a 'complete' inner product (scalar product). 

We have the fluxes ]u, 3 pi , and —of with the corresponding forces 

= ^^(^) (28) 

=TV(^^^) (29) 
L = Vv. (30) 



In the last step of the derivation, the thermodynamic fluxes have to be determined. Following On- 
sager, we first make the basic ansatz of linear non-equilibrium thermodynamics, i.e. we postulate 
a linear dependence of the fluxes on all the forces (if symmetry permits), i.e. 

j/3 = Ei^/^a-f" (31) 

a 

with L^^ = Jjp^ (Onsager relation, resulting from microscopic reversibility). Specifically, we get 
for the fluxes 

^)+rL„„.v(i^-''- 

i)+TL„,..v(«^-''^ 
= -LggiVv (34) 

Note that there is no linear coupling between the momentum flux and the thermodynamic forces 
corresponding to temperature and concentration gradients. However, when discussing the total 
energy for systems with large gradients in the density pi we included quadratic terms in the density 
gradient. For consistency, a nonlinear term, quadratic in the forces fp^, should be added to relation 
(l34l) resulting in 

- -L„ : Vv - Q^^, : (v ^^^) (v i^^) . (35) 
The additional term is related to irreversible aspects of the dynamics of the diffuse interface and 



TL„„ ■ V ( - ) + TL„p^ ■ V ( ^^^^^ 1 (32) 
= ^Lp..- V ( ^ ] + TL,,,, ■ V ( ) (33) 



can be seen as a generalization of the term proposed in the conclusion of Ref. 
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that reads 



S. = /3(VC)(VC), (36) 



where P is an undetermined empirical coefficient and is the difference between the chemical 
potentials of components 1 and 2 for an inhomogeneous equilibrium which in our terms is C = 
/id — ^^Pi as is discussed below, after Eq. (l57l) and also in Ref. T^. Q^^ corresponds to a te nsor 
of interfacial viscosities. Related issues are discussed for sharp interface theories in Ref. 



102. 



The Lq,^ are tensors of various orders: L„„ is of order 2, whereas Lgg and are of order 4 



—SPi 

34 



(i.e. in the general case there are 3 = 81 components: viscosities). Assumption of an isotropic 
medium significantly reduces the number of coefficients.— 

Considering small interfacial viscosities only, we neglect the corresponding terms and finally 
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get 



3 Pi 



^Api - /id 
T 

^Api - /id 
T 



-CKV ■ v) - r/ 



Vv + (Vv 



ll(V ■ v) 



(37) 
(38) 
(39) 



where ( and r/ are the dynamic and shear bulk viscosity, respectively, and the kinetic coefficients 
have the units [A;i] = kgm/s'^, [r/] = [Q = [k2] = kg/ms and [fcs] = kgs/m'^. The governing 
equations are now obtained by introducing the fluxes (l37l) to (|39l) into the corresponding transport 
equations. Introducing (l38l) into (fTOl) we get for the transport of pi 



^Pi + V ■ (piv) + V 



-^VT + A;3rV 



|^_^Api - /id 



T 



and for momentum transport feeding (l39l) and (l27l) into ([8]) we have 



P- + PV.VV 
-a(V-v)-r/ 



V- U(Vpi)(Vpi) + U - ePiAp 







^(Vpi)^ 



(40) 



Vv + (Vv)T - -I(V ■ v) 



(41) 



Assuming that /i2 and /id are given, to obtain a closed set of equations we still need an equation 
for the evolution of the temperature field. We obtain it by multiplying the transport equation for 
entropy (l24l) with T, and ex press ing the entropy via the Gibbs relation in a similar way as Batchelor 



(p.35ff andp. 136ff of Ref. 



IQQ) 



After reordering Eq. (|24)) yields 



-TV 



+ R. 



(42) 



All terms on the r.h.s. are already known. We next determine the l.h.s. terms. Consistently with 
Gibbs relation (fTSi) one can express the entropy as a function of T, V, Ni and Using the second 
law of thermodynamics we write 



dQ = TdS = CvdT + T — 



dV + T 



dS 



dNi + T 



T,V,N2 



dS 



dNo 



dN2, (43) 



2 / T,V,Ni 



where Q stands for the heat supplied to the system. Fixing A^^i and N2 the last two terms vanish. 
Expressing volume V in terms of T, p, Ni, and N2, assuming fixed A^i and N2 yields 



dV 



dV 



dT 



— dT+\ — 



dV 



dp 



dp 



T,Ni 
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Substituting into Eq. (1431) we obtain (assuming fixed A^i and A^2) 
TdS 

dp. 



dT + T(—) (—) dp 



(44) 



Introducing the thermal expansion coefficient (3 = y {^)pNNi comparing the respective 
prefactors of dT in the two lines of Eq. ((44)) we rewrite Eq. (|43l) as 

dS 



dQ = TdS = CydT + '^'^ ^S^ dV + T 



(3V 



dNi 



dNi+T 



dm 



(45) 



T,V,N2 / T,V,Ni 

Assuming that the material constants Cp, Cy and (3 are known we next focus on the last two terms 
of (|45l) . Consistently with Eq. (fTSl ). we write for the Helmholtz free energy, F = U — TS, of a 
local fluid element V in thermodynamic equilibrium 



dF = d{U - TS) = -SdT - pdV + fiidNi + fi2dN2. 
Partial differentiation of (|46|) with respect to A^i and T gives 



(46) 



dS 
dm 



dfii 
dT 



T,V,N2 



V,N, 



d 


(dF 


dm 




-i 


' dF 


dT \ 


.dm 



and 



V,N. 



T,V,N2 



(47) 



Ty,N2 / y 



N,, 



respectively. Identifying the mixed 2nd derivatives we obtain the so called Maxwell relation for S 

and fii 



dS 

dm / T,V,N2 



dT 



Analogously, we obtain the Maxwell relation for S and /t2. Using them we rewrite (l45l) as 



TdS = CvdT + ^^^dV-T(^ 

dT 



dm-T 



dT 



dm.. 



(49) 



\ "^-^ / VM \ "^-^ / V,Ni 

Note that there exist other ways to express the heat. The one chosen here is advantageous because 
the equation explicitly contains dV. This easily allows to consider the incompressible case by set- 
ting dV = (see below Section HTCl) . The dependencies of the chemical potentials on temperature 
will be also given. 
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Next, let us rewrite local Gibbs and Gibbs-Duhem relations, Eqs. (l20l) and (|2T1) . in terms of the 
density of the Helmholtz free energy, /(p, pi, T, t) = u — Ts 



df = —sdT + fi2dp + Pddpi, 
P= -f + l^2p + Pdpi- 



(50) 

(51) 



using the above introduced procedure for obtaining local Gibbs (l20l) and Gibbs-Duhem (|2TI) rela- 
tions, we derive from Eq. (|49l) , the relations for the volume densities of the extensive quantities 

'dp2\ , rj. f dpd 



Tds = cydT — T 

Cp - Cv 



dT 



dp — T 



p,pi 



Ts 



/5 



T 



dp2 
dT 



p-T 



dT 

dpd 
dT 



dpi 



and 



(52) 
(53) 



p,pi ^ ^ p,pi 

respectively. As above, p2 = ^2^0/^2 and pd = PiNa/Mi — p2Na/M2 are related to densities 
not to particle number. Note that we also changed the notation for partial derivative with respect 
to T regarding the chemical potentials p2, Pd as defined by (|50l) . 

Dividing finally Eq. (|52l) by (it, identifying d/dt with the material derivative D/Dt, substitut- 
ing into Eq. (I42l) . using Eqs. ([3]) and (flOl) and reordering we obtain the transport equation for the 
temperature field 

DT Cr,-cv— „ fdpd 



T 



dT 



V-jp, = -TV 



P:Pl 



i?. (54) 



C. Model-H - bulk equations 

Next, we simplify the coupled equations for temperature, momentum and volume density of 
component 1 by assuming a fluid with constant density p (which implies incompressibility V ■ v = 
0), in an isothermal setting (constant T). Further, we express the density pi in terms of a mass 
concentration ci = pi/p and obtain from Eq. (|40l) the convective Cahn-Hilliard equation 

dci 



dt 



+ V • Vci + V • [MiV (a,, Aci - Pd) 



0, 



(55) 



where we introduced pd = PPd, = P^^ and Mi = k^/p'^. Note that ci is dimensionless and 



[fjcj = mkg/s^ [pd] = kg/ms^ and [Mi] = sm^kg 
The momentum equation (|4T]) reduces to 



p-^ + pv • Vv = -V ■ |(Tci(Vci)(Vci) + p - (TciCiAci 2 



^(Vci)2] ij+r^Av, (56) 
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where the mechanical pressure is given by the local Gibbs-Duhem relation (|2T]) . i.e. p = —u + 
IJL2P + Ts + jidCi- We emphasize that p as the mechanical pressure for a homogeneous material in 
the thermodynamic equilibrium is a locally defined variable that should not depend on any gradient 
or derivative. However, for simplicity we introduce an effective pressure p^s^ that incorporates all 
terms in the square brackets in equation (|56l) . Equations (|55] ) and (|56l ) are normally called model- 

In the literature model-H is presented in various forms. Especially, the momentum equation 
is written in different ways. Most differences arise from different definitions of the pressure PeSi 
introduced in the last paragraph. In the following we indicate how to translate the different formu- 
lations and point out 'irreducible' differences. 

The review by Anderson et al.^^ gives as transport equation for the momentum in a binary 
mixture (their Eq. (17b) with (16a) and (19)) in our notation) 



P 



-v + vVv 



-V- 



(PA-^(Vci)^)I + (T,,(Vci)(Vci) +r/Av 



(57) 



However, their equation (20) for their chemical potential /i^ indicates that they do not follow their 
Eq. (10), but already absorbed additional terms into their Our formulation coincides with theirs 
identifying their — o"ci(Vci)^/2 and our pcfTi- The difference in the formulation arises because 
Anderson et al. use /i^ in place of our /id — cTc^ Aci = fi in the thermodynamic pressure definition, 
i.e. their /ic = /i is the chemical potential for an inhomogeneous equilibrium. Then also their Eq. 
(21) corresponds to our (|55l) . 

Jasnov and Vinals^ present two forms for the momentum equation 



P 



-v + vVv 



(58) 



[their Eq. (2)] and 

I d ] ^ 

p —V + V ■ Vv = -Vpjv + ?7Av - ciV/i. (59) 
\_dt J 

[their Eq. (2) with the replacement described in the last paragraph of their appendix]. They also 
use /2 instead of /id- We introduce different symbols pjv and pjv for the respective pressures. 
The second form can be obtained from ours taking into account V • [pcffil + crci(Vci)(Vci)] = 
ci V/2 + Vp. The first form just follows from integration by parts and redefining the pressure again: 
Pjv = Pjv - ci/i. 

The form of model-H presented in the review by Hohenberg and Halperin^^ gives a momen- 
tum equation (their Eq. (5.1b)) that agrees on the first view with the second form of Jasnov and 
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Vinals (|59| ). However, they dropped the pressure term, i.e. in the limit of constant concentration 
their model does not reduce to the Navier-Stokes equations. 

Finally, we rewrite model-H in terms of the difference of concentrations c — Ci — C2 — 2cx — 1. 
Introducing new parameters ac = cFci/^ and M = 4Mi and specifying the chemical potential 
/id = 2dcf{c), where /(c) is the concentration dependent part of the local free energy, results in 

dtc + v-Vc = -V • {MV [(TeAc - dj{c)]}. (60) 

and 

dv 

+ pv ■ Vv = -V ■ {a,(Vc)(Vc) +pefrl} + r?Av (61) 

where 

Pes = P - CX,{C + 1) Ac - y(Vc)2. (g2) 

Fixing /(c) to be a symmetric double well potential, Eq. (l60l) corresponds to the convective Cahn- 
Hilliard equation studied, for instance, in Ref.-^^^*^. The energy will be further discussed in 
SectionlYl Because of its importance for the boundary conditions (see Section [nil) we also give 
the stress tensor 

r = -Peffl - ae(Vc)(Vc) + r/ (Vv + (Vv)^) , (63) 

where {p — Pcs)l — o-c(Vc) (Vc) represents the Korteweg stress.-^^ The pressure Pcs can be cal- 
culated from the Poisson equation 

Apeff = -fXe(VV) : [(Vc)(Vc)] - p(Vv) : (Vv). (64) 

In the literature the various formulations of model-H are mainly used to describe the behavior 
of bulk flows . "^"^'"^^1"^^'^° Systems confined between rigid plates are considered in some cases^^*^^ 
assuming (i) the diffuse interface is far away from the plates, and (ii) the walls are neutral with 
respect to the two components. However, the role of energetically biased plates and the evolution 
of a free surface of the binary mixture have to be understood in their interaction with the bulk flow 
to be able to describe an evolving free surface film on a solid support. The necessary boundary 
conditions are discussed next. 
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III. BOUNDARY CONDITIONS 



A. Concentration 

For the concentration field the boundary conditions were discussed in connection with a purely 
diffusive transport for a system confined by rigid plates. ■^^iZiti^ Assuming the velocity is zero at the 
rigid substrate (no-slip condition, see Section ITlIB I) the conditions for the full model-H are similar. 
We have zero diffusive flux through the substrate {z = 0) 

d,[a,Ac~dJ{c)] = (65) 

and obtain in the general case an evolution equation for the concentration (see Appendix lAl) 

dtC + Y-Vc = -M-[-aAc - a-Aiic + dj-{c)] (66) 

where Ay = Vy ■ Vy and V|| = {d^, dy). Here, however, we will focus on a surface energy that 
(i) does not depend on concentration gradients (cr^ = 0) and (ii) relaxes instantaneously to its 
equilibrium value (A/^ oo). 

At the free surface {z = h{x,y,t)] one has the condition of zero diffusive flux through the 
moving surface, i.e. n ■ j^^ = with j^^ as defined in Eqs. (flOl) and (l38l) . 

is the normal vector of the free surface. The change from the total flux j^^ [Eq. ©J to jp^ = 
— piv exactly accounts for the transformation into the frame moving locally with the surface. 
One gets 

n- V(a,Ac-aj(c)) = 0. (68) 

The second condition is in the general case again an evolution equation for the concentration field 
on the boundary as derived in Appendix |Al The evolution equation is valid in the local comoving 
frame, i.e. 

dtC + v-Vc = -M+[a,{n ■ V)c - a+A^c + dj+{c)], (69) 

where = ■ Vs and the surface nabla operator is defined as = (I — nn) ■ V. In the 
following we assume as above (t+ = and M+ — oo. We will drop the respective terms after 
the non-dimensionalization in Section |IVl Note, however, that the units of the surface parameters 
differ from the ones of the corresponding bulk parameters: [M"^] = s/kg, [af] = kgm^/s^, 
= N/m = kg/s2. 
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B. Velocity 



The boundary conditions for the velocity fields are the no-slip and no-penetration condition at 
the solid substrate {z = 0) 

V = 0, (70) 
and the force equilibrium at the free surface (z = h) 

(r-T,,,) -n = -7(c) nV-n + V.7(c) (71) 

Note that V ■ n corresponds to the curvature of the free surface. We assume that the ambient air 
does not transmit any force (r^jr = 0)- The first term on the r.h.s. of Eq. (TtTI) corresponds to the 
Laplace or curvature pressure whereas the second one represents a Marangoni force tangential to 
the interface and results from the variation of the surface tension along the surface caused normally 
by a solutal or thermal Marangoni effect. As shown in the Appendix |A] these terms can be derived 
from a minimization procedure. 

For a pure Navier-Stokes problem the Marangoni term is often modeled as a linear dependence 
of the surface tension on concentration or temperature. Here, however, one has to use a condition in 
accordance with the interface energies introduced when discussing the boundary conditions for the 
concentration field, i.e. at Eq. (l69l) . For cr+ = and A/+ oo [see Eq. (l69l) 1 the surface tension 
7(c) corresponds to /'^(c) plus a constant (reference tension 70, see below SectionlV]). For (t+ 7^ 
the surface tension depends as well on concentration gradients 7 = 7(0, (Vsc)^), a concept that 
has not yet been followed in the literature. Considering a finite Af + would correspond to a 7(0, t), 
i.e. to a dynamical surface tension characterized by a relaxation time towards its equilibrium value. 
Both complications will not be considered further in the present paper. 

The boundary condition (ItTI) is of vectorial character, i.e. three scalar conditions are derived by 
projecting it onto n, ti, and t2, respectively, where 

^ (l,0,g./i) ^ (0,1,^^^ 

are the (non orthogonal) tangent vectors. The resulting scalar conditions 

— (Tc (n ■ Vc)^ — Peflf + 2r7 n • (Vv) • n = —7(c) V ■ n 
~ac (tl ■ Vc)(n ■ Vc) + r^ti ■ (Vv + (Vv)^) ■ n = ti 
-ac (ta ■ Vc)(n • Vc) + r/ts • (Vv + (Vv)^) ■ n = ta 
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(72) 
(73) 

■ V7(c) (74) 

■ V7(c) (75) 



correspond to equilibria of normal and tangential forces, respectively. 

At the free surface one has furthermore the kinematic condition, i.e. the prescription that the 
surface follows the flow field 

dth = n-v ^1 + (V||/i)2 (76) 
which can be written in a more compact form as 

(dth) • n = V • n (77) 

where vector h = h{x,y,t) is tracking the free surface. 



IV. NON-DIMENSIONALIZATION 



Next, we non-dimensionalize the bulk equations and the boundary equations in 2 steps: (i) 
introduction of abstract scales for velocity, pressure, length, concentration and energy density that 
leads us to a set of dimensionless numbers; (ii) introduction of problem specific length and velocity 
scales to obtain a minimal set of dimensionless numbers valid for the problems without external 
driving studied here. 



A. Abstract scales 



Introducing scales 



dimensionless 


Scale 


dimensional 


t' 


T„ = l/U 


t = Tyt' 


x' 


I 




v' 


u 




p' 


p 


p = Pp' 


c' 


c 


c = Cc' 


fV) 


E 


f{c) - Ef\d 









one obtains after dropping the primes the dimensionless bulk equations 

9tc + V • Vc = -Ts V • {V [KoAc - dj{c)]} . (79) 
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and 



Ps 
where 



- V • {Ko' ( Vc) ( Vc) + peff 1} + ^ Av (80) 

K.e 



Peff = P - Ko'(c + 1) Ac - iKo'(Vc)2. (81) 
We defined the dimensionless numbers 

Reynolds number Re = 

V 

Korteweg number 1 Ko = — r- — 



Korteweg number 2 Ko' 



PE 

(y. 



ME 

Time scale ratio Ts = 
Pressure scale ratio Ps = — ^ (82) 

We propose the name 'Korteweg number' because both of them are related to the Korteweg 
stresses. The Korteweg numbers can be seen as 'bulk Marangoni numbers'. For the determination 
of the energy scale E see SectionlVl 

The scaled boundary conditions for the concentration field at both interfaces are the no-flux 
condition 

= n - V[KoAc-9J(c)]. (83) 
and the evolution equations for the concentration at the surface 

dtC + w-Vc = -Ts^ [Ko n ■ Vc - Ko^ A,c + En^9c/^(c)]. (84) 

For the substrate one sets n = (0, 0, —1). 

At the free surface the conditions for the normal and tangential forces are 

Ps 

-Ko' (n ■ Vcf - peff + 2— n • (Vv) ■ n = -S 7(c) V ■ n (85) 

Ps 

-Ko' (ti ■ Vc)(n ■ Vc) + — ti ■ (Vv + (Vv)^) • n = S ti ■ V7(c) (86) 

Re 

Ps 

-Ko' (t2 ■ Vc)(n ■ Vc) + — t2 ■ (Vv + (Vv)^) • n = S ta ■ V7(c) (87) 

Re 

respectively, where 7 is the dimensionless surface tension referred below in Section |V] as 7'. The 
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dimensionless numbers are either given above or listed next 

Boundary Korteweg number Ko — 



PE 

, EfM^ 
Boundary time scale ratio Ts = 

U (_y 

Boundary energy number En = - — 

IE 

Surface tension number ^ ~ Jp ^^^^ 



B. Specific scales 

For relaxational settings, i.e. systems without external driving forces one might specify scales 
based on the 'internal' diffusive or convective transport. Assuming very viscous liquids and taking 
into account that all structure formation will be driven by the decomposition process, it is conve- 
nient to base all scales on the diffusive processes only. Fixing Ko = 1 and Ts = 1 length and 
velocity scales become 

E^ ^ = TC 

respectively. Choosing a pressure scale based on the energy density scale 




l=J^C and U = ^^, (89) 



P = E (90) 

identifies the two Korteweg numbers, i.e. Ko'=Ko= 1. The specific forms of Reynolds and Pres- 
sure numbers are then 



MEp , pM^E pM^E'^ 



respectively, i.e. Ps/Re= r]M/l'^C^ = rfME / acC^ ■ The dimensionless numbers related to the 
boundaries become 

Note that Ts"*^ oo for — > oo and Ko"*^ = for cr^ = (see discussion in Section UlI)) . 



, E^ E^ 7o 7o 
En^ = — = —TTo and S = = (93) 

IE al'^E^^C IE al'^E^^C 

The relation of S and the 'classical' Marangoni number is discussed in the next section. 



20 



V. LOCAL ENERGIES 



For the local bulk and surface energies we use a simple polynomial, an approximation valid 
near the critical point. However, it is straightforward to introduce other expressions derived using 
Flory-Huggins or more advanced theories.- For the bulk energy a symmetric quartic potential is 
used 

m = /„ - + (94) 



corresponding to the nondimensional form (Eqs. ([78]) ') 

nC) = \ - 1)' + const (95) 

with E = hC^ and C = ^Ja/h. For the surface energies of the two interfaces we use the respective 
quadratic expressions 

/^(c) = 70^= + + yc2. (96) 

Note that in the framework of model-H for a film of binary mixture the surface energies /~(c) 
and /^(c) correspond to the concentration dependent surface tensions of the liquid-solid and the 
liquid-gas interface, respectively. This implies that f'^{c) is responsible for a linear (6+ = 0) or 
nonlinear (6+ 7^ 0) Marangoni effect. The surface energies and 7^^ are the respective reference 
surface tensions at c = 0. 

Using the reference surface tension of the free surface as a scale for both interfaces, i.e. = 
E~ = 'Jq = 7o, we arrive at the nondimensional expressions 

f'Hc) = 4 + + (97) 
7o 2 

and identify 7' = /'"""(c'). The dimensionless parameters 

= —— and = —— (98) 

7o 7o 

describe preferential adsorption of one of the species at the interface and changes in the interaction 
between the species at the respective interfaces. Inspecting Eq. (|86l) it becomes clear that the 
'classical' Marangoni number for a linear Marangoni effect is Ma = a+S. The corr espondin g 



105 



number for a quadratic Marangoni effect is Ma2 = (compare, for instance, Ref. 
Furthermore, we can now specify 

En+ = En" = ^ = S. (99) 

LE 
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This implies that the boundary conditions for concentration and momentum transport at the free 
surface are intrinsically coupled. Note, finally that these considerations only apply for Ts^ — »• 
oo and Ko^ = 0. See discussion in Sections Hill andHVl In the following we only work with 
dimensionless quantities and drop all primes. 

We are now equipped with a complete model to investigate a wide variety of systems involving 
decomposing mixtures with free surfaces. Although, the boundary conditions in Section Hill are 
written for a film on a solid substrate they can be easily adapted for free standing films, i.e. for a 
film with two free surfaces. Also droplets of a mixture on a solid substrate can be studied if the 
given model is supplemented by a condition at the contact line such as a concentration dependent 
equilibrium contact angle. This will be the scope of future work. 

To understand the evolution of the surface and concentration profiles of a decomposing film we 
next analyse (i) the homogeneous and vertically stratified base state solutions, (ii) the transversal 
instability of the base state solutions that lead to the experimentally observed film profiles and 
concentration patterns, and (iii) the full non-linear time evolution. Part (i) will be studied in the 
remainder of the present paper, part (ii) forms the content of the accompanying paper, and part 
(iii) will be presented in a planned sequel. 

VI. BASE STATES 

The understanding of the behaviour of a thin film of a mixture on a solid substrate has to be 
based on an analysis of the base state solutions. For a film on a horizontal substrate without further 
driving forces parallel to the substrate the base states are quiescent, i.e. the velocity of the fluid 
mixture is zero. We distinguish two types of quiescent base states: (a) completely homogeneous 
flat film and (b) horizontally (transversally) homogeneous but vertically stratified film. 

A. Completely homogeneous film 

A completely homogeneous film of arbitrary thickness h{x, y) = ho, with arbitrary concentra- 
tion c(x) = Co and with quiescent fluid vq = corresponds to a base state solution of the system 
(r79l)-(l97l) only if there exists no energetic bias at the solid-liquid interface or the free surface, i.e. 
without any linear or nonlinear Marangoni effect: = 6^ = 0. The corresponding effective 
pressure is Pcs = 
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For energetically biased interfaces the boundary conditions for the concentration field are only 
fulfilled if dcf^ico) = 0, i.e. for cq = —a'^/b'^ = —a^/b^. For finite a+, a~ this is from the 
experimental point of view a very unlikely case. Here we will not pursue it further. However, for 
a"*" = a~ = a homogeneous film of a symmetric mixture (cq = 0) represents a base state for 
any fe+ and b^ . This case corresponds to a purely quadratic Marangoni effect. Experimentally, it 
is not a very common case but was studied in hydrodynamics for films of alcohol solutions^ and 
also as a problem of purely diffusive demixing in a gap. The latter case was anal ysed in detail in 



Refs. 



74 



75 



and will serve as a benchmark for our linear stability analysis in Ref. 1 107 



B. Vertically stratified, horizontally homogeneous film 

Depositing a thin film of a mixture on a solid substrate it is to expect that processes that lead 
to a vertical stratification are much faster than processes that lead to a horizontal structuring if the 
film thickness is similar or below the length scale of bulk decomposition. The vertically stratified 
films may on a larger time scale undergo a further horizontal structuring. The finally emerging 
horizontal length scales and structures can be understood from the 'short-time' vertical layering. 
Therefore we focus next on a systematic investigation of steady layered films. 

A flat layer (h = h^) of a quiescent fluid mixture (vq = 0) represents a base state if the 
vertical concentration profile c = Co{z) is a steady solution of the classical one-dimensional non- 
convective Cahn-Hilliard equation 

dtc = d,,[d,,c-dj{c)]. (100) 

and the boundary conditions (i) = dz [dzzCo ~ f^c/(co)] {at z = and z = h) and (ii) = 
[±9zCo + S9c/^(co)] ('— ' at z = and '+' dX z = h). Taking into account (i) one has to solve 
the bulk equation dzzCQ — <9c/(co) + /^i = with boundary conditions (ii). The constant of 
integration Ki represents the dimensionless chemical potential for an inhomogeneous equilibrium 
as discussed after (Eq.lSTi). 

The remaining equations and boundary conditions are fulfilled withpcfr = Pca{z) = —{dzCoY + 
const, i.e. the layers are completely characterized by co{z). In the following we determine families 
of solutions in terms of concentration profiles for (i) energetically neutral or non-biased surfaces, 
(ii) symmetrically biased surfaces, (iii) antisymmetrically biased surfaces, and (iv) asymmetrically 
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FIG. 1: Branches of steady vertical concentration profiles for energetically non-biased (neutral) surfaces 
(o^ = 6^ = 0) and a symmetric mixture (c = 0) in dependence on the film thickness h. Shown are (a) the 
L2-norm and (b) the energy E of co{z). Selected corresponding profiles are given in Fig.|2] 5=1, and the 
symbols are explained in the main text. 



biased surfaces. Thereby we characterize the concentration profiles by their energy 



E 



r + r + / [{d.cY + /(c)] dz 

'o 



and the L2-norm 



ll^cll 



h Jo 



\c(z) — c]'^ dz 



(101) 



(102) 



where c is the mean concentration. Note that E should only be used to compare films of ide ntical 

iotI 



h and c. The profiles are determined using numerical continuation techniques detailed in Ref. 
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1. Energetically neutral surfaces 




c 



FIG. 2: Selected concentration profiles corresponding to Fig. [T] Note that each of the solutions has a 
'twin' obtained by c ^ — c that also corresponds to an allowed profile. This degeneracy may be lifted by 
energetically biased surfaces (depending on the symmetry z ^ — z, see below). The n = 1/2 profiles 
corresponding to the energy minimum for the respective film thickness are shown as heavy (red online) 
lines. 

The results for the trivial case of energetically neutral solid substrate and free film surface are 
given in Fig.[T]for a symmetric mixture, i.e. the case of zero mean concentration c = 0. Shown are 
the L2-norm and the energy E per film area. Fig. [2]presents selected concentration profiles. The 
base states for a film correspond to selected solutions of the one-dimensional bulk Cahn-Hilliard 
equation. A multiple of the period has to be equal to the film thickness. At both, the substrate and 
the free surface, the profile has a minimum or a maximum. This allows to classify the obtained 
solution branches by the number of periods n. The simplest stratified films correspond to half a 
period (n = 1/2), one period {n = 1), one and a half period (n = 3/2) and so on. 

Note that the solutions with an integer n are symmetric with respect to a reflection at the plane 
z = ho/2, i.e. co{z) = co(/io — z). We call them in the following 'z-reflection-symmetric'. 
They are accompanied by a twin solution with identical L2-norm and energy obtained by an in- 
version of concentration: cq{z) —cq{z). On the contrary, the solutions with a non-integer n 
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are antisymmetric with respect to a reflection at the plane z = /io/2, i.e. co{z) = — co(/io — z) 
('2; -reflection-antisymmetric'). The resulting second solution has naturally identical L2-norm and 
energy, and can also be obtained by an inversion of concentration. We will also call it the twin 
solution. 

The bifurcation diagram is not very involved. For all film thicknesses there exists the trivial 
homogeneous solution with | |5c| | = = (broken line in Fig. [T]). Non-trivial branches bifurcate 
successively at hi = in/kc where i = 1,2,3. . . and kc = a/— i9cc/(co) = a/1 — Scg corre- 
sponds to the critical wavenumber for the linear instability of the homogeneous solution c = Cq 
of Eq. (|100l) . For the symmetric mixture considered in Fig. [Hone finds hi = in. Furthermore, 
all characteristics like Ei[h] or \ \Sc\ \.i[h] of all branches n = i with i > 1 can be mapped onto the 
characteristics of the n = 1/2 branch. For example, for the energy one has Ei[h] = Ei/2[h/{2i)]. 
Note, however, that the bifurcations are degenerate because as discussed above two twin solutions 
related by symmetry bifurcate at once. 

A thin film in an experiment will tend towards the constellation with the minimal energy (see 
Figs. \T\h and |2l), i.e. for h < n the homogeneous layer and for h > n the stratified layer with 
n = 1/2. The multilayer constellations with n > 1 may, however, appear as transients as they are 
saddle fixed points in phase space that attract time evolutions from a certain basin of attraction and 
repel them consecutively into the few unstable directions (for a more extensive discussion of that 

Q 

concept in connection with dewetting on heterogeneous substrates see Ref. I49|) . 

2. Symmetrically biased surfaces 

The presented rather detailed description of the steady states for energetically neutral surfaces 
will help us to understand the involved behaviour for biased surfaces. Allowing for arbitrary 
linear (a^, a+) and quadratic (6^, 6+) energetic biases opens a four dimensional parameter space 
additionally to the parameter 'film thickness'. We give an overview of the system behaviour by 
focusing on a linear bias (b^ = 6+ = 0), and by using several special ratios a^/a~. In this way 
we obtain a 2d parameter space spanned by a+ and h. 

In the present section we assume that the two surfaces energetically prefer the same compo- 
nent with equal strength (a+ = a^), i.e. we have symmetrically biased surfaces. Figs. |3]and|4] 
show characteristics of solution branches in dependence of the bias for fixed film thickness and 
in dependence of the film thickness for fixed bias, respectively. Corresponding solutions between 
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1 



n=l/2 




FIG. 3: Branches of steady vertical concentration profiles for linearly symmetrically biased surfaces {a~^ = 
a~ = 0) and a symmetric mixture (c = 0) in dependence of bias a for film thicknesses as given in the 
legend. Shown are (a) the L2-norm and (b) the energy E of cq{z). Selected corresponding profiles are given 
in Fig. |5] The branch numbers n for the different film thicknesses are given using different fonts: /i = 3.5 
bold, h = 5 italic, h = 10 normal. Lines and labels are of corresponding colors (online). S = 1, and the 
symbols are explained in the main text. 
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Figs. [Hand [3] are marked by filled symbols in the L2-norm plots. Hollow symbols indicate cor- 
respondences between Figs. [3] and HI Concentration profiles for a+ = 0.2 corresponding to the 
hollow symbols are given in Fig.[5](a) whereas panel (b) gives profiles for a large bias of a+ = 0.6. 
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FIG. 4: Branches of steady vertical concentration profiles in dependence of film thicknesses for linearly 
symmetrically biased surfaces (a+ = = 0.2, 6^ = 0) and a symmetric mixture (c = 0). Shown are (a) 
the L2-norm and (b) the energy E of co{z). Selected corresponding profiles are given in Fig. [51 S = 1, and 
the symbols are explained in the main text. 



First, we focus on Fig. [3l We introduce branch names indicating the 'non-biased branch' 
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FIG. 5: Selected concentration profiles corresponding to Figs. [3] and |4] for (a) a+ = 0.2 and (b) a+ = 0.6 
sorted by branch number as indicated in Fig. |4l Film thicknesses are indicated in the legends. The profiles 
corresponding to the energy minimum for the respective film thickness and bias are shown as heavy (red 
online) lines. The minuscule 'b' in the n = 1 panel for o"*" = 0.6 denotes the profile on the side branch of 
the n = 1 branch (cp. Fig. (3]). 
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(Fig. [T]) they are emerging from. This convention does not correspond to actual properties of 
the concentration profile. For instance, the profiles on the n = branch in Fig. [3] are not homo- 
geneous any more. For small film thicknesses n < h < 2n only two solutions exist at a+ = 
corresponding to two branches for increasing a+ > 0. Thereby the n = [n = 1/2] branch is 
unstable [stable]. Further increasing the bias the two branches approach each other. For h = 3.5 
the n = 1/2 branch terminates in a supercritical bifurcation on the n = branch. For h = 5 the 
stable n = 1/2 branch first undergoes a saddle-node bifurcation turning unstable before it finally 
terminates in a subcritical bifurcation on the n = branch. Beyond the bifurcation the n = 
branch is stable in both cases. 

For i7c < h < {i + l)7c one finds i solutions at a+ = exemplified in Fig. [3]for /i = 10, 
where 4 solutions exist. Increasing a+ > one finds, however, 5 emerging branches because the 
degeneracy of the n = 1 solution at a+ = is lifted by the energetic bias. See the discussion of 
symmetries above in Section IVIB 1[ 

The 2; -reflection- symmetric solutions (integer n) have at a'*" = a 'twin' -solution obtained by 
c{z) — > —c{z) that reacts differently when imposing a symmetric energetic bias, i.e. the bias lifts 
the degeneracy and two distinct branches are generated like, for instance, in Fig. |3]for /i = 10 
and n = 1 (see also profiles in Fig. [5]). Note that the n = branch is a special case with out 
degeneracy at a+ = (trivial solution at a+ = 0). On the contrary, the 'twin' -solutions of the 
^-reflection-antisymmetric solutions (non-integer n) do not react in a different way to a symmetric 
bias, i.e. their degeneracy is not lifted. 

Coming back to the case h = 10 we see that when further increasing a+ most branches end 
in bifurcations. One branch finally survives for large bias corresponding to a layer of 'liquid +' 
enclosed by two layers of 'liquid — ' that is preferred by both surfaces for a+ > 0. This implies 
that depending on the strength of bias the energetic minimum corresponds to qualitatively different 
stratifications - bilayer (n = 1/2) structure for small a+ and a sandwich trilayer (n = 1) structure 
for large a+ (see heavy (red online) lines in Fig. [51 

Note finally that the diagram is symmetric w.r.t. a+ — a+. Focusing on the branch of lowest 
energy that represents the solutions selected by the system we see that the role is taken for small 
[large] a+ by the n = 1/2 [n = 1 or n = 0] branch. It is intuitively clear that a strong symmetric 
bias will suppress the 2;-reflection-antisymmetric solutions. The alternative view of fixing a+ and 
changing h is given for a+ = 0.2 in Fig. [Hallowing for a better comparison with Fig. [TJ From this 
representation it becomes clear that for a+ > the 2 branches emerging from the n = 1 solution 
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-n=l/2 




FIG. 6: Branches of steady vertical concentration profiles for linearly antisymmetrically biased surfaces 
{a~^ = —a~,b^ = 0) and a symmetric mixture (c = 0) in dependence of bias a for film thicknesses as 
given in the legend. Shown are (a) the L2-norm and (b) the energy E of cq{z). Selected corresponding 
profiles are given in Fig.|7J 5 = 1. The branch numbers n for the different film thicknesses are given in 
different fonts: h = 3.5 bold, h = 5 italic, h = 10 normal. Lines and labels are of corresponding colors 
(online). 
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'break off the n = branch a.th = 27r. A similar process occurs at all h = 2i7r for integer i. 
3. Antisymmetrically biased surfaces 

In contrast to the preceding section, here we assume a"*" = —a^, i.e. the two surfaces ener- 
getically prefer different components. The preference is, however, equally strong. We focus on 
a+ > 0, i.e. the free surface prefers the c < component. The case a+ < is related by symmetry. 
Figs. [6] and [8] show solution branches in dependence of the bias and of film thickness, respectively 
(in analogy to Figs.[3]and |4l). Selected corresponding profiles are given in Fig.|71 



n=0 n=l/2 n=l n=3/2 




FIG. 7: Selected concentration profiles corresponding to Fig.|6]for a"*" = 0.2 sorted by branch number as 
indicated in Fig. [6l Film thicknesses are indicated in the legend. The profiles corresponding to the energy 
minimum for the respective film thickness and bias are shown as heavy (red online) lines. 

In contrast to the case of symmetrically biased surfaces we find that for small a+ > two 
branches emerge from the n = 1/2 and n = 3/2 solutions but only one from the n = 1 solution. 
Here, the degeneracy of the solutions at a+ = is only lifted for the z-reflection-antisymmetric 
solutions (non-integer n), but not for the 2 -reflection- symmetric solutions (integer n). In the for- 
mer case one of the twin solutions at a+ = is favored by the antisymmetric bias whereas the 
other one is disfavored, i.e. they decrease and increase their energy with a+, respectively (Fig.[6]b). 
One of the n = 1/2 solutions is the only one that 'survives' for large bias a~^. It is furthermore this 
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solution that corresponds to the energy minimum for all a"*" (see heavy (red online) lines in Fig. IT]). 

The alternative view of fixing a+ and changing h is given for a"*" = 0.2 in Fig. [8l Contrary 
to section IVI B 21 for > the n = branch is 'broken off at h = 2mx by the respective 2 
branches emerging from the non-integer n solutions. 
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FIG. 8: Branches of steady vertical concentration profiles for linearly antisymmetrically biased surfaces 
{a^ = —a^ = 0.2, = 0) and a symmetric mixture (c = 0) in dependence of film thickness. Shown are 
(a) the L2-norm and (b) the energy E of co{z). Selected corresponding profiles are given in Fig.|7] S = I. 
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4. Asymmetrically biased surfaces 
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FIG. 9: Branches of steady vertical concentration profiles for linearly asymmetrically biased surfaces {a^ > 
0, a~ = 6^ = 0) and a symmetric mixture (c = 0) in dependence of bias a+ for film thicknesses as given 
in the legend. Shown are (a) the L2-norm and (b) the energy E of cq{^z). 5* = 1. The branch numbers n for 
the different film thicknesses are given in different fonts: h = 3.5 bold, h = 5 italic, h = 10 normal. Lines 
and labels are of corresponding colors (online). 



As an intermediate case compared to the two preceding sections, we focus next on a = 
and a+ > 0, i.e. the substrate is energetically neutral whereas the free surface prefers the c < 
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component. Note that the cases a~ = and a+ < 0, a~ > and a"*" = 0, a~ < and a"*" = are 
related by symmetry. 

Figs. |9] and [To] show solution branches in dependence of the bias for fixed film thickness and in 
dependence of the film thickness for fixed bias, respectively. Here the bias lifts all degeneracies 
existing for a+ = 0, i.e. from each solution at a+ = emerge two branches (beside the n = 
branch). Correspondingly, Fig. [10] shows that for a+ > the 2 branches emerging from every 
integer and non-integer n solution 'break off' the n = branch at all h = 2mi. 

The branch of lowest energy is here for all a+ and h the n = 1/2 branch, i.e. a simple two 
layer structure. This is, however, by no means a general result but depends on the specific asym- 
metry chosen. For a strong bias that is only slightly asymmetric (like, for instance, = a+ + A 
with A ^ a+ the branch of lowest energy will still be the sandwich structure discussed in Sec- 
tion JvHU 

C. Non-symmetric mixtures 

We have seen that the case of a symmetric mixture c = can well be used to understand the 
basic solution structure for stratified films. However, it has to be kept in mind that a symmetric 
mixture represents a rather special case. Experimental systems will normally consist of non- 
symmetric mixtures with c ^ 0. In the present section we give selected results for the general 
case. 

For small c ^ the branch structure for non-biased surfaces is given in Fig. [TT] for c = 0.25. 
It is qualitatively equivalent to the one for a symmetric mixture (cp. Fig. [T]). Note, however, that 
in our normalization the energy for the trivial homogeneous film now depends linearly on film 
thickness because the bulk energy density for c is not zero any more. For larger c the primary 
bifurcations become subcritical. Eventually the trivial solution becomes linearly stable. It is, 
however, metastable, because finite perturbations may trigger a nonlinear instability. For those c 
no primary bifurcations exist. All branches of stratified solutions continue towards infinite thick- 
ness. The bifurcation diagrams for biased surfaces become quite involved for the subcritical and 
metastable case and will be discussed elsewhere. 

Focusing on the case of supercritical primary bifurcations shown in Fig. \TT\ next we discuss 
the influence of symmetrically biased surfaces presented in Fig. [121 The general form of the 
bifurcation diagram for a+ > is qualitatively very similar to the case of a symmetric mixture 
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FIG. 10: Branches of steady vertical concentration profiles for linearly asymmetrically biased surfaces 
(a+ = 0.2, = 0, = 0) and a symmetric mixture (c = 0) in dependence of film thicknesses. Shown 
are (a) the L2-norm and (b) the energy E of co{z). S = I. 

(Fig. |3]). However, the symmetry w.r.t. a+ — a+ does not hold anymore. It is replaced by a 
symmetry w.r.t. (c, a+) (— c, — a+). 

Also for antisymmetrically biased surfaces one finds similar bifurcation diagrams for non- 
symmetric (Fig. [13]) and symmetric (Fig. (6]) mixtures. The antisymmetry of the boundary con- 
ditions implies that the symmetry w.r.t. a+ — a+ also holds for the non-symmetric mixture. 
Finally, in the asymmetrically biased case shown in Fig. [14] all degeneracies at a+ = are bro- 
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FIG. 1 1 : Branches of steady vertical concentration profiles for non-biased surfaces (a^ = 6^ = 0) and 
non-symmetric mixtures (c = 0.25) in dependence of film thickness h. Shown are (a) the L2-nomi and (b) 
the energy E of cq{z). 5 = 1. 

ken as in the case of a symmetric mixture (Fig. [9l). Furthermore, now also the symmetry w.r.t. 
a+ — a+ is broken. 

VII. CONCLUSION 

We have proposed a dynamical model that describes the coupled decomposition and profile 
evolution of a free surface film of a binary mixture, a process frequently encountered in coating 
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FIG. 12: Branches of steady vertical concentration profiles for linearly symmetrically biased surfaces (a"^ = 
a^,b^ = 0) and a non-symmetric mixture (c = 0.25) in dependence of bias a+ for film thicknesses as given 
in the legend. Shown are (a) the L2-norm and (b) the energy E of co{z). S = 1. The branch numbers n for 
the different film thicknesses are given in different fonts: h = 5 italic, h = W normal. Lines and labels are 
of corresponding colors (online). 

and structuring processes. An example is a thin film of a polymer blend on a solid substrate un- 
dergoing simultaneous phase separation and dewetting. We have based our approach on model-H 
coupling transport of the mass of one component (convective Cahn-Hilliard equation) and mo- 
mentum (Navier-Stokes-Korteweg equations). We have used the framework of phenomenological 
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FIG. 13: Branches of steady vertical concentration profiles for linearly antisymmetrically biased surfaces 
(a"'" = —a~, = 0) and a non-symmetric mixture (c = 0.25) in dependence of bias for film thicknesses 
as given in the legend. Shown are (a) the L2-norm and (b) the energy E of co{z). 5 = 1. The branch 
numbers n for the different film thicknesses are given in different fonts: h = 5 italic, h = W normal. Lines 
and labels are of corresponding colors (online). 

non-equilibrium thermodynamics to derive a generalized model-H coupling transport equations 
for momentum, density and entropy in the framework of phenomenological non-equilibrium ther- 
modynamics. Then we have discussed the individual contributions representing, for instance, an 
internal Soret effect and interface viscosity. The model has been simplified for isothermal setting. 
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FIG. 14: Branches of steady vertical concentration profiles for linearly asymmetrically biased surfaces 
(a+ 7^ 0, a~ = 0, = 0) and a non-symmetric mixture (c = 0.25) in dependence of bias a+ for film 
thicknesses as given in the legend. Shown are (a) the L2-norm and (b) the energy E of co{z). S = 1. The 
branch numbers n for the different film thicknesses are given in different fonts: h = 5 italic, h = W normal. 
Lines and labels are of corresponding colors (online). 

vanishing interface viscosity and internal energies resulting from a setting close to the critical point 
of demixing. A comparison with literature results has clarified the issue of defining pressure and 
chemical potential. 

To facilitate the description of a free surface profile we have introduced boundary conditions 
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at the solid substrate and the free interface. It has been explained that the incorporation of hy- 
drodynamic flow even in the case of extremely slow creeping flow is a necessary precondition for 
the description of evolving surface deflections. It has been shown that the dimensionless num- 
bers entering the boundary conditions for the Cahn-Hilliard and the Korteweg-Navier-Stokes are 
closely related. Therefore they can not by any means be chosen independently of each other. 
After non-dimensionalization we have analysed possible steady base state solutions for laterally 
homogeneous films of decomposing mixtures. In doing so we have distinguished vertically ho- 
mogeneous and vertically stratified films. It has been shown that a plethora of stratified solutions 
exist that can be mapped, ordered and understood using continuation techniques and symmetry ar- 
guments. The obtained systematics will form the basis for the analysis of the lateral stability of the 
base states undertaken in an accompanying paper. '"^ In an Appendix lAl we have used variational 
calculus to independently confirm the boundary conditions for the static limiting case. 

Our results on vertical stratifications can be compared to a broad variety of experimental data 
on static film structures. However, they can also be used to interpret transitions observed in slow 
time evolutions. Most of the results on vertical layering reviewed in^ can be explained at least qual- 
itatively. Most static vertical layerings observed in thin films of polymer blends either correspond 
to two-layer or sandwich-like three-layer structure that we have found to be the only structures of 
lowest energy depending on the energetic bias of the surfaces. 

The interesting case of the evolution of a relatively thick (500 nm) decomposing d-PMMA/S AN 
blend film^ is presented in Fig. 16 of Ref. 4. The vertical profile develops from a homogeneous 
film to a two-layer structure, with d-PMMA collecting at the solid substrate. However, the path 
to equilibrium passes through a sandwich-like three-layer structure. This can be easily understood 
from the solution structure presented for antisymmetrically biased surfaces in Fig. |6l There the 
n = 1 (sandwich-like three-layer structure) has a lower energy than the n = Corn = 3/2 
solution but a higher energy than the final n = 1/2 solution. The n = 1 solution corresponds to a 
saddle in phase space, i.e. it attracts time evolutions of a broad range of initial conditions and then 
expels evolutions in its only unstable direction directing the evolution towards the final two-layer 
structure. Which 'saddle solutions' will be involved in a time evolution depends on the wavelength 
of the fastest linear mode. For a similar discussion for dewetting films on a heterogeneous substrate 



see Ref. 
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section 3. 



In another experiment using a PEP/d-PEP blend film^ presented in Fig. 2 1 of Ref. Q it is shown 
that the equilibrium layer structure changes from d-PEP/PEP/d-PEP to d-PEP/PEP by varying the 



41 



substrate surface energy. This corresponds in our idealized setting to a transition from symmetri- 
cally biased (or unbiased) surfaces to anti- or asymmetrically biased surfaces. The here observed 
change in the layering of lowest energy corresponds well to the experiments. 

To conclude, the present work has presented a complete model to investigate a wide variety 
of systems involving the coupled evolution of surface and concentration profiles of free surface 
films of a decomposing mixture. It can be used to analyse vertically stratified but horizontally 
homogeneous films and their evolution in time. This includes layer inversions of two-layer sys- 
tems with diffuse or sharp interfaces that can not be described by two-layer models for immiscible 
liquids. ■^2iiLd24^ The dynamics of such an inversion is quite complex. For instance, for a sub- 
strate/liquid 1 /liquid 2/air two-layer structure it involves the transient formation of drops of liquid 1 
at the liquid 2/air interface. The typical distance of those droplets can now be calculated using a 
linear transversal stability analysis of the unstable two-layer profile. For antisymmetrically biased 
surfaces it exists, for instance, for /i = 5 up to a+ ~ 0.3 (Fig. (6]). 

Our model also allows to study the emergence of permanent lateral structures using a transver- 
sal stability analysis of the stratified layers^^ or a simulation in time. Especially, it allows to 
quantify the influence of hydrodynamic transport on the dynamics. 

Note that the model can be adapted for several related problems. Although, the boundary 
conditions in Section [ni] are posed for a film on a solid substrate they can easily be adapted for 
free standing films, i.e. for a film with two free surfaces. Also drops of a mixture on a solid 
substrate can be studied if the presented model is supplemented by a condition at the contact line. 
The latter will be the scope of future work. Also the used model of local bulk and interface energies 
can easily be replaced by more realistic functions as discussed in Ref. 4l. Choosing parameters that 
correspond to a stable mixture also the dynamics of mixing in a free surface film can be studied. 



APPENDIX A: VARIATIONAL APPROACH 



This appendix uses variational calculus to derive the static limit of the bulk equations and 
boundary conditions from the underlying Helmholtz free energy functional of Cahn-Hilliard'- 
type for a two-dimensional film of binary mixture. The three dimensional case will be presented 
elsewhere for a more general setting. The free energy functional 

F[c{x, z), h{x)] = F^[c{x, z), h{x)] + F,[c{x, h{x)), h{x)]. (Al) 
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consists of a bulk part Fb and a surface part Fg defined as 

-co Jo 



dzdx — Ad 



ci dA - CiA 



A 



dA-A 



h{c)ds. 



(A2) 
(A3) 



dQ 



The second and third integral in Fb are taken over the same area as the first one, the symbol Q is 
introduced for convenience. The symbol dQ denotes the boundary of the domain of integration 
Q. We assume that the surface free energy does not depend on (Vc)^, i.e., the bulk free energy 
gradient term is not amended at the free surface. Such a contribution and its physical consequences 
will be discussed elsewhere. The Lagrange multipliers Ad and A enforce mean concentration of 
component 1, ci = (c + l)/2, and total area of the domain to have the prescribed values Ci and 
A, respectively. Later, we will discuss their relation to the local chemical potentials /i2, Atd and 
mechanical pressure p. 

To vary F with respect to all possible degrees of freedom in a transparent way we define small 
changes of the functions c and h as 



h{x-^ a) = h{x) + Oi({x), 
c{x, z; a) — c(x; a) — c(x) + Q;r^(x), 



(A4) 
(A5) 



where C and t] are arbitrary admissible functions and x = (x, z). The variation of h and c cor- 
responds to the operation da\a=o, i-e., 5h{x) — ({x), 5c(x) = r]{x). However, not only the 
local concentration and the interface position are varied. Varying the latter also implies that fluid 
elements have to vary their position due to convective motion. 

The variation of the bulk contribution to the free energy functional can be written as 



SF. 



oo rh{x) 
oo JQ 



dzdx — 5 Ad 



cidA - CiA 



SX 



dA-A 



a 



+ 



^(Vcf + /(c)-AdCi-A ax)dx, 



(A6) 



where we used 5Vc = V5c. Integrating by parts the term containing V77, we get 



5Fb 



00 ph{x) 



-00 -'0 



-ae(Ac)r; + ( - ^ ) 



dzdx — 5Ad 



cidA - CiA 



-SX 



dA-A 



n 



+ ® aciyc) ■ nrjds + 

'dQ 



-(Vc)2 + /(c)-AdCi-A C{x)dx. 



(A7) 
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Next, we turn our attention to the surface contribution Fg. concentrating on the most interesting 
top part of the boundary df2, i.e. the free surface. We denote the corresponding part of Fg as F^°p 
and write it as 

F.*°P= / fJcU.)]^dx, (A8) 



where 



ds 

dx 



ds 

/s[c(Xs)]— dx, 



(A9) 



For the point Xs(a;) = (x, h{x)) at the free surface and its variation we have 

Xs(x; a) = (x, h{x) + a({x)), 
(5x3(x) = (0,C(x)). 



(AlO) 



Then 5c(xs(x)) = Vc(xs) ■5xs+?7(xs). Using eqs. (IA9I) and Xs(x) we have expressed the integrand 
of (IA8I) as a function of x. The variation of F*°p is 



dcfs ^ 5c(xs) + /s(c)5^ }dx 



j |9Js(Vc . + r^)^ + Uc)d,ht ■ Vcjrfx, (All) 



where we used Eq. (IA9I) and applied 

9s^(x) rfC(3;) c^^; rfC(3;) 



t-vc- 



(A12) 



(ds/dx) dx ds ds 

The last step is correct if fields are only defined at the surface (here (, but valid also for h, n, t, 
etc.) are interpreted as being defined everywhere with values independent of z. 

Next, integration by parts has to be applied to eliminate derivatives of the variations. One uses 

ja{t-Wh)ds = J a{t-Vb)^dx = - J bV ■ t^a dx = - J b {t ■ Va) ds. (A13) 

Assuming laterally periodic or localized structures, surface terms resulting from the integration by 
parts are zero here. We obtain 



(5Ff P = / <^ dJ^iVc ■ (5x, + r/) - t ■ V 
Performing derivatives and substituting from (lAlOl ) we arrive at 



C>ds. 



(A14) 



5Ff P = jsjs [{d.c)C + [(t ■ V/s) d^h - ^C^ds, (A15) 
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where 

dccxh 

- -JiTWW' ^^^^^ 

is the curvature. Note that k is positive for a convex surface of the fluid. Terms with {dcfs)C can 
be simplified as follows 



(9,c)— - (t-Vc)9./i 



-C=(9./.)(n.Vc)-C (A17) 



resulting in the final expression 

5Ff P = { [(n • Vc) + /s «:] + (^c/s)^?} (is. (A18) 

The variation of the contribution of the free energy at the bottom part of the boundary df2, 
denoted by 5F^°^, can be obtained as a special case of 5F^°p. In (|A11)) we consider Xs(x) = 
(x, 0), = 0, ds/dx = 1 resulting in 

5F':°'= / 9J,r/da:. (A19) 



Next, writing 5F = 5Fy, + + 5F^°^ = 0, one is in principle ready to extract governing 

equations and natural boundary conditions of the problem. Inspecting the form of (IA7I) . (|A18I) . and 
(|A19I) one notes that we obtained two scalar boundary conditions on the free surface as the pref- 
actors of arbitrary admissible functions i], in the boundary integral of the stationarity condition 

6F = 

(Ten ■ Vc + dj, = 0, (A20) 
^(Vc)^ + /(c) - AdCi - A + (n ■ Vc) dj, + U = 0, (A21) 

respectively. In order to obtain force boundary conditions, we need to express our variations (, 77 
in terms of the virtual displacements because mechanical forces are energetically conjugated to 
them. 

For this purpose, we introduce the variation of the position of a fluid element due to convective 
motion. The varied Euler coordinates x of a fluid element specified by its material (Lagrange) 
coordinates X can be expressed as 

x(X;a) = X + ax(X), (A22) 
5x(X) = x(X). (A23) 
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where x = iXx, Xz) is an arbitrary admissible displacement vector. Both coordinate systems 
coincide for a = 0. The free surface has to follow the fluid as is expressed by the kinematic 
condition Eq. dTT] ). This introduces a dependency between the variations C and x, i.e. 



dx 

6h{x) = C{x) = -d^h{x)xx{x) + Xz{x), or C ■ n = C ^ = X • n, (A24) 

Next, we consider a fluid element identified by its referential position X. The varied concentration 
at this element can be expressed using (|A22|) as 

c(x; a) = c(X + ax(X)) + ar]{X + ax(X)), (A25) 



being consistent with (IA5I) . Consider for a moment that no diffusion is active. In that case the 
concentration c(x; a) changes only due to convection described by X- I^i consequence, the con- 
centration at the arbitrary but fixed fluid element X should remain constant, i.e. 

^ {c[X + axm + ar][X + ax(X)]} = Vc[X+ax(X)]-x(X) + r/[X+ax(X)] = 0. (A26) 

We denote the variation r] that satisfies this condition as r]co (convective) and the remaining part as 
?7di (diffusive) variation. In consequence, we have 

11 = 1100 + i]di = - Vc -X + ildi (A27) 

where ?7di is a variation independent of x because it is caused by a different physical process. 

Using Eqs. (IA27I) and (|A24I) . we write the stationarity conditions for the variation 5F = 5F^ + 
SFg. Prefactors of the variations in the bulk and at the free and bottom surface give the Euler- 
Lagrange equations of the problem. The prefactor of r^di in the bulk integral gives 

-a,Ac + dJ-^ = 0, (A28) 



i.e. the correct static limit of the 2d version of Eq. (1601) . From this equation we deduce that 
— 2(TcAc + 2dcf = — cTci Aci + /id = Ad is the chemical potential for a heterogeneous equilibrium 
discussed in Section IllCi The prefactor of x gives the same equation as (IA28I) . The surface 
integrals yield as the prefactor of rj^i 

CTell- Vc + = 0, (A29) 

i.e., the static limit of the 2d version of Eq. (I66l) with a~ = and of Eq. (|69l ) with (t+ = 0. On the 
free surface the vectorial prefactor of x 

(n ■ Vc) (9 Js) n - ( Vc) dj, + f^^n - a,n ■ {Vc) ( Vc) 



+ 



Ly(Vc)^ + /(c)-AdCi-A 



n = (A30) 
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gives, using I — nn = tt and reordering, 

-a,n ■ (Vc)(Vc) - pefrn = tt ■ V/s - fsKn (A31) 

with 

Pes = {dcf - aeAc)(c + 1) - y (Vc)2 - /(c) + A, (A32) 

where we used Eq. (IA28I) for Ad We proceed to identify the quantities in peS- Using (ISTI) . we 
substitute the free energy density /(c) in (IA32I) and simplify to 

Pes=p- I^2P - oAc + l)Ac - y (Vc)=' + A, (A33) 

taking into account that 2{dcf)ci = jldCi = ^dPi- The pressure pcs in (|A33I) coincides with the 
one defined in (|62l) in case that 

A = /i2p. (A34) 

We conclude that the Eq. (IA31I) gives the tangential and the normal force equilibrium conditions 
at the free surface and corresponds to the static limit of the boundary conditions (TtTI) with (|63l) . 
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